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High transmission plateaus exist widely in the logarithmic transmission spectra of localized sys- 
tems. Their physical origins are short chains of coupled-localized-states embedded inside the lo- 
calized system, which are dubbed as "short necklace states". In this work, we define the essential 
quantities and then, based on these quantities, we investigate the short necklace states' properties 
statistically and quantitatively. Two different approaches are utilized and the results from them 
agree with each other very well. In the first approach, the typical plateau-width and the typical 
order of short necklace states are obtained from the correlation function of logarithmic transmission. 
In the second approach, we investigate statistical distributions of the peak/plateau-width measured 
in logarithmic transmission spectra. A novel distribution is found, which can be exactly fitted by the 
summation of two Gaussian distributions. These two distributions are the results of sharp peaks of 
localized states and the high plateaus of short necklace states. The center of the second distribution 
also tells us the typical plateau-width of short necklace states. With increasing the system length, 
the scaling property of typical plateau-width is very special since it almost does not decrease. The 
methods and the quantities defined in this work can be widely used on Anderson localization studies. 

PACS numbers: 42.25.Dd, 72.15.Rn, 78.67.Pt, 42.25.Bs 
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I. INTRODUCTION 

The most extraordinary phenomenon of wave trans- 
port in random media is Anderson localization 1 , which is 
shared by various systems, such as, electronic, photonic 
and accoustic systems 2 -. In Anderson localized system 
the ensemble average of logarithmic transmission (InT), 
not transmission (T), is additive with system length 
l£&. Naturally, the localization length can be defined as 
£ = —2L/ (InT). In strongly localized systems(L ;§> £), 
the transmission is generally small, but also shows large 
fluctuations^. Those extreme large T values will domi- 
nate the transmission of localized systems. In the past 
years, the study of the physical origins of those large 
T values yielded abundant results for understanding the 
transport phenomena in random systems and the Ander- 
son localization 5 - - — . 

Previously, it was pointed out by M. A. Azbel et al* 5 - 
that the resonant transport through those localized states 
lying near the system center can contribute very large T 
values, being of order unity. Although such a mecha- 
nism can generate high resonant transmission peak 5 -^, 
its contribution to the transmission is vanishingly small 
in strongly localized system because the resonant peak is 
exponentially sharp (~ e~ L ^)^&&. Laterly, Pendrj* 7 - and 
Tartakovskii 8 independently predicted a kind of quasi- 
extended state, called the necklace state(NS). The NS is 
formed through the coupling between nearly-degenerated 
localized states which are evenly distributed in the sys- 
tem. Despite the spatial overlaps between localized states 
are small, because of the degenerate coupling, the NS con- 
tributes very wide transmission "mini-bands" , which can 
dominate the transmission of strongly localized systems. 
The NSs have been demonstrated experimentally in pho- 



tonic system s 10 ! 11 . Their fundamental characters and 
statistical consequences to the transmission are widely 
studie d 12 ! 14 " — . Recently, it is demonstrated that the NSs 
have evident contribution to the short-time transport of 
wave package 1 ^— and the dynamics of fluctuations of lo- 
calized waves22. More profoundly, the number of NS can 
increase dramatically as approaching the Anderson tran- 
sition point, which strongly supports a modes-coupling 
induced quantum percolation scenario for the Anderson 
localization-delocalization transitio n 24 ! 25 . 

Even though NSs can contribute very large transmis- 
sion, their formation has rigorous conditions, i.e., the 
nearly-degenerated localized states are required to be 
evenly distributed inside the systero L 14 ' 16 ! 17 . In gen- 
eral the localized states with similar frequency are not 
such well distributed in a specific random configuration. 
In such cases, ideal NS a 14 ^ 15 crossing the whole config- 
uration are not formed. But eventually, those nearly- 
degenerated states which are also spatially close to each 
other couple together, forming short chains of coupled- 
localized-states embedded inside the configuration. Since 
those coupled-localized-states have similar properties as 
the NS, we call them short necklace state (SNS) 26 . The 
SNS also manifests as high transmission plateau which 
has significant transmission contribution. In transmis- 
sion spectra the SNS is hardly to be distinguished from 
isolated localized states since the valley between coupled 
peaks seems very low. But it can be clearly identified 
from the logarithmic transmission spectra, for instance, 
see the coupled-peaks in Fig.l. Because of the strong 
coupling between the neighboring localized states inside 
the SNS, on top of the plateau there are sharp peaks and 
valleys, which correspond to un-smooth changes of the 
transmission phase between the coupled-peaksi 4 -. 



In contrast to the NSs, the SNSs have special 
properties 2 ^ that (1) the SNSs exist widely in every ran- 
dom configuration while the NS only occurs once in mil- 
lions of random configurations with large length(i 3> 
£); (2) the plateau- width of the SNSs only depends on 
the coupling strength between the neighboring localized 
states and is insensitive to the system length. With 
these properties, SNSs are superior qualitatively com- 
pared with the NSs. However, how to quantitatively 
characterize statistical properties and scaling behaviors 
of SNSs are still unexplored topic. Even more, some es- 
sential quantities of SNS study, such as the "half-width" 
of SNS plateaus in InT spectra, need to be defined since 
there is no obvious physical quantity in previous studies 
can directly describe SNS. 

In this paper, we study the statistical manifestation 
of the SNS using physical quantities which are defined 
in the logarithmic transmission spectra. To characterize 
SNS, we find out two different approaches, whose results 
can be compared with each other. The first approach is 
based on the correlation functions Ci n T (defined in the 
InT spectra) and C't (defined in the transmission coef- 
ficient t spectra). From C;„t, we show that there is a 
typical frequency correlation range of logarithmic trans- 
mission, which is explained as the typical SNS plateau- 
width. From this frequency range and compared with 
the results of Ct, we can find the most-probable order 
of SNS. The second approach is from the direct mea- 
surement of peak/plateau- width in InT spectra. We find 
the statistical distribution of peak/plateau- width is quite 
abnormal and can be fitted very well by the summation 
of two Gaussian distribution functions, where the pri- 
mary Gaussian center gives the typical width of resonant 
peaks of localized states while the secondary one gives 
exactly the typical width of the SNS plateaus. Excellent 
agreements are found between two approaches. The de- 
pendencies of SNS's properties on the scaling parameter 
L/£ are also studied. We find the plateau-width of SNS 
almost does not depend on L/£ while the peak-width of 
localized states decays exponentially with L/£, indicating 
the SNSs have more significant transmission contribution 
in longer systems. The essential quantities defined in this 
work also provide new ways for further quantitative study 
on statistical properties of the transmission of Anderson 
localized systems. 

The rest of this paper is organized as follows. In sec- 
tion II, we introduce our model and basic properties of 
the InT spectra. In section III, we present our numer- 
ical results and theoretical analysis of three correlation 
functions: the correlation of transmission Ct, the corre- 
lation of logarithmic transmission Ci n x and the field cor- 
relation Ct- We show that the half-width of Ci n T gives 
the typical plateau-width of SNS. The most-probable or- 
der of SNS can be obtained by comparing Ci u t and Ct- 
In section IV, we study statistical distributions of the 
peak/plateau- width (defined in the InT spectra). The 
probability distributions are calculated from very large 
numbers of samples and can be fitted very well by the 
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FIG. 1. (Color online) Typical logarithmic transmission spec- 
tra (a) and corresponding transmission spectra (b) of L — 10£ 
systems. Red bold line: £ = 34.13a (W = 0.6); black thin line: 
£ = 215.6a (W = 0.24). The SNS marked in (a) is enlarged 
in (b) for clear observation. 



summation of two Gaussian distributions, where the first 
Gaussian center gives the typical width of resonant peaks 
of localized states and the second one gives the typical 
width of the SNS plateau. In section V, we discuss the 
system length dependency of the SNS. Finally, a sum- 
mary of this work is given in section VI. 



II. 



THE MODEL 



We study ID random stacks composed of binary- 
dielectric layers(called A and B) with thicknesses d,A — 
d,B = 500nm and with refractive indices n^ = 1.0 
and us = 3.0 * (1 + W-y), where 7 is a random num- 
ber uniformly distributed in [—0.5,0.5] and W gives 
the randomness strength. Such a periodic on average 
model is an optical counterpart of the Anderson model 
of electronic systems and is widely used for localization 
studie a 12 ' 14 i 15 ' 19 i 25 . The transmission coefficients of op- 
tical waves are calculated by standard transfer matrix 
method 25 . Without randomness, W — 0, the system ex- 
hibits a pass band in frequency range (89.13, \2A.lTHz). 
Our study will focus on the range (100.0, 103.2T.ffz), 
where the localization length is almost a constant for a 
certain randomness W(0.2 < W < 0.6) 25 . The localiza- 
tion length is calculated by £ = —2L/ (InT), where (InT) 
is averaged from the 10 6 configurations satisfying L>( 
(L«500£)- 

Fig. 1(a) shows two typical logarithmic transmission 
spectra with the same L/£ = 10 but with different 
localization lengths, £ = 215.6a(black thin) and £ = 



34.14a(red bold). We can see they have similar aver- 
age height, being approximately — 2L/£ = —20. After 
averaged over a large number of configurations, the InT 
spectrum becomes a flat line exactly falls on the mean 
value (InT) = —20. This is a natural result according 
to the single parameter scaling theorem that InT fol- 
lows the Gaussian distribution with mean value — 2L/£. 
An interesting phenomenon in the InT spectra is that 
there are some clear plateau-structures, such as the one 
marked by arrows in Fig. 1(a). From the wave inten- 
sity distributions one can find those plateaus are formed 
by the SNSs(i.e., from the degenerate coupling between 
some spatially close localized states) embedded inside the 
configuration 26 . Since the peaks of localized states in 
the spectra are extremely sharp (~ e~ L ^), the fluctua- 
tion of InT is dominated by those plateaus. More im- 
portantly, these SNS may be essential for understanding 
the Anderson transition phenomenon 26 . After observed 
a large number of InT spectra for fixed £ and L/£, we 
find the plateaus appear with some intuitionistic regu- 
larities, (1) most plateaus have similar frequency width 
for a configuration; (2) the number of peaks on each 
plateau are similar. For example, for most plateaus of 
the £ = 34.14a system shown in Fig.l, the plateau width 
is about 0.3 ~ OATHz and the number of peaks are 
about 3^4. In the following we will focus on the sta- 
tistical properties of these SNSs and try to characterize 
them quantitatively. 

To quantitatively study the SNS, we need to find 
proper physical quantities. Fig.l shows that the plateaus 
of SNSs can hardly be distinguished from localized states 
in a T spectrum (see Fig. 1(b)), but can be clearly identi- 
fied by the plateaus in the InT spectra (Fig. 1(a)). So we 
try to define physical quantities in the InT spectra for the 
SNS study. In the following, we will define and study the 
correlation function of InT spectra in our first approach. 
Then, we will define the InT peak/plateau width-T and 
study the statistical distribution of T in our second ap- 
proach. 



III. 



CORRELATION FUNCTIONS 



According to the standard definition of mathematics, 
we define the correlation functions of T and InT in fre- 
quency domain as: 



Ct(Aw) 



Cj nT (Aw) 
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where T is the transmission coefficient, er is the stan- 
dard deviation a(x) — \/((x — (a;)) 2 ) and Cov is the co- 
variance Cov(x,y) = (x — (x)) (y — (y)). Physically Ct 
implicates the frequency correlation of transmitted wave 
intensity. Its Fourier transformation corresponds to the 



dynamical response at the outgoing interface, which have 
been intensively studied in weakly localized systems 27 . 
If two frequencies are on the same transmission peak, 
Ct is close to unity. Otherwise Ct will close to zero. 
Hence the half-width of Ct is usually used to character- 
ize the linewidth of localized states^ 7 -. Similarly, on the 
InT spectra, since InT of two frequencies on the same 
plateau are much larger than the mean value (InT) and 
contribute a significantly large Ci n T- We expect Ci u t is 
close to unity when Au is smaller than a typical plateau- 
width and decays to zero when Aw is larger than a typical 
plateau- width. 

Meanwhile, we also study the field correlation function, 



C t (Aw) 



< T u > + < Toj+Aw > 



(3) 



where t is complex transmission coefficient and T = tt* . 
The complex amplitude of incident wave is chosen to 
be Ei n (x = 0) = 1 so that at the outgoing interface 
t = E(x = L) = \E\e 1 ^, where E is the complex elec- 
tronic field and <j> is its phase. Similar to Ct, the phys- 
ical implication of Ct is simply the frequency correla- 
tion of transmitted field, i.e., the field correlation (Ct is 
the intensity correlation). The definition of Ct is the 
same as that in reb 7 -, which could be negative valued, 
depending on the averaged phase difference A(j> between 
w and w + Aw. Generally, when Aw crosses a resonant 
peak, the phase of t jumps tt— — . Hence the phase dif- 
ference between t{ui) and t(u> + Aw) is approximately 
Atp = cf>(u))-4>(u)+Au) « 7r. Then £ w £* +Aw + £* Waw = 
l-E^H-E'id+Aidl ■ 2cos(A(j)) becomes negative and so that C t 
is negative. When Aw gets to the value which typically 
contains two peaks (plus one average distance between 
two peaks), Ct will be positive since A<fi w 2-7T, and so 
forth. Hence when increasing Aw, we expect Ct(Aw) 
to oscillate between positive and negative. Each time 
Ct(Aw) changes its sign, the average number of resonant 
peaks in Aw increases one. So Ct(Aw) provides a way 
to check the average resonant peak number in a certain 
frequency range. 

In our calculations we set the original frequency point 
w = lQQTHz without lose of generality. The correlation 
functions calculated from a very large number (2 x 10 7 ) of 
configurations for different £ are shown in Fig. 2(a), where 
the ratio L/£ is fixed at 10. The solid, dashed and dot- 
ted curves respectively represents Ct , Ct and Ci n T ■ The 
red curves marked by crosses correspond to the smaller £ 
system and the black curves marked by solid dots corre- 
spond to the larger £ system. In Fig. 2(b) the transverse 
axis is shown in logarithmic scale. 

Let us first discuss Ct- Fig. 2(a) shows Ct is a very 
singular function at the origin Aw — > 0. The linewidth 
Aw of localized states, corresponding to CV(Aw) = 0.5, is 
larger in the smaller £ systems. This is in consistent with 
the observation in Fig.l that the linewidth of resonant 
peaks in smaller £ system is larger. Detailed data gives 
that the halfwidth of the larger £ system (black solid 
curve) is about 0.382GHz and the smaller £ system (red 
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FIG. 2. (Color online) (a): Correlation functions of L = 10£ 
systems calculated from 2 x 10 7 configurations. Red: £ = 
34.13a; Black: £ = 215.6a. For each color, the correlation 
functions from the top down are respectively Ci n T, Ct and 
Ct- (b): Same as (a) but with Alo shown in logarithmic scale, 
(c): Detailed view of the oscillations of Ct(£ = 34.13a). The 
half- width of Ci„t, fis, gives the typical plateau- width of 
SNS. The first minimum of Ct, £2i, gives the typical width of 
a resonant peak. 



solid curve) is 2A83GHz. This is also quantitatively in 
consistent with the observation from the InT spectra. 

Ci n T is a much smoother function than Ct in the 
Acj — > limit. Take the £ = 34.13a system (red curves) 
for example. Ct rapidly falls to zero near ln(Au>) m 24. 
Ci n T exhibits a plateau at the origin, then drops to 0.5 
at In(Aoj) « 26.6, i.e., its half-width is about OMITHz, 
being much larger than that of Ct- Such contrast reflects 
the different geometry properties between the T and InT 
spectrum. On the T spectrum, Ct falls to zero typically 
when two frequencies are not on the same peak. Since 
the peaks of localized states are exponentially sharp, Ct 
decreases rapidly as increasing Aw. However, on the InT 
spectrum, there exist many plateaus formed by SNSs, 
as shown in Fig. 1(a). Those plateaus are usually higher 
than the average transmission background (InT). When 
two frequencies are on the same plateau but not the same 
peak, it will still contribute a significantly large value to 
CinT- Hence the halfwidth of Ci u t, which is denoted as 
Qg i- 71 this paper, approximately gives the typical plateau- 
width of SNS in InT spectrum. 

With the help of Ct , we can find more detailed informa- 
tion of SNS, such as the average SNS order, which is the 
average number of coupled localized states in one SNS. 
Similar to Ct, Ct also shows sharp singularity at the 
origin and falls quickly as increasing Aw. Interestingly, 
Ct shows some oscillations around Ct = in the region 
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FIG. 3. (Color online) Probability distribution of the 
peak/plateau-width F(see the text for definition) measured 
from logarithmic transmission spectrum. The distribution 
function can be fitted very well by the summation(red solid 
curve) of two Gaussian functions (dashed curves). The pri- 
mary(left, blue) Gaussian center corresponds to the typi- 
cal width of resonant peaks of localized states. The sec- 
ondary (right, green) Gaussian center corresponds to the typ- 
ical plateau- width of SNS. 



where Ct falls to nearly zero, see Fig. 2(c). As discussed 
earlier in this section, those oscillations can be under- 
stood from the 7r-phase jumps of resonant peaks of the 
localized states. The different order minimal/maximal 
points of the Ct correspond to the frequency ranges which 
can accommodate certain number of resonant peaks. We 
denote the Aw at the first minimum of Ct as fix, which is 
the the typical frequency width of a resonant peak. We 
also denote nth order minimal/maximal points VL n as the 
average frequency range which can accommodate n reso- 
nant peaks. So Ct is the ruler of the number of resonant 
peaks in a frequency range. With VL n in our mind, we 
can see that the mean width of SNS plateaus, obtained 
by CinT, can accommodate about 3^4 resonant peaks, 
as indicated by the green arrows in the Fig. 2(c). Actually, 
this average number of resonant peaks in a SNS agrees 
very well with our direct observation of many spectra. 



IV. PEAK/PLATEAU- WIDTH STATISTICS 

To make a test and verify of the characters of SNS 
obtained from Ci n T and Ct, we next try our second ap- 
proach, i.e. direct measure of the peaks/plateaus- width 
in a large number of InT spectra. Since the ensem- 
ble average of InT spectra- (InT) is well defined, it is 
natural to define the peaks /plateaus-width T as the fre- 
quency interval where InT is always higher than (InT). 
From direct observation of the InT spectra, one can find 
that (InT) crosses lots of sharp peaks of localized states 
and fewer SNS plateaus. More precise statistical results 
should be obtained from a large number of realizations. 
We find that the statistical distribution of T is extremely 
skewed. The reason is that in Anderson localized sys- 
tems the peak-width of localized states is exponentially 
small and the width of coupled-peaks also scales exponen- 
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FIG. 4. (Color online) (a): Ct (lower four blue curves) and 
Ci n T (upper four red curves) for different L/£ systems. For 
each set of four curves, L/£ = 6, 10, 14, 18 from left to right. 
The inset gives detailed view at the halfwidth of Ci n T- (b): 
The typical peak- width of localized states Hi (marked by blue 
arrows in (a)) as a function of L/£. (c): The typical plateau- 
width of SNS ils (marked by red arrows in (a)) as a function 
of L/£. Hi decays exponentially as increasing L/£ while ils 
is insensitive to L/£. 



tially with the system length£ii2£. This is very similar to 
the probability distribution of the dimensionless conduc- 
tance g, which is extremely skewed (nearly log- normal). 
Early study^ on Anderson localization has shown that 
it is better to study the additive quantity-Zng, which is 
nearly Gaussian distributed. Similarly, instead of T, we 
will study the InT distribution, which is likely Gaussian 
distributed. 

We have measured T in 10 spectra with high frequency 
precision that can distinguish each resonant peak in the 
L = 10£ systems^,. The probability distributions of InT 
for two different £ systems are shown in Fig. 3. Take the 
£ = 34.13a system as the example. (All the following 
discussions also apply to the £ = 215.6a system.) The 
probability distribution is roughly Gaussian-like, where 
a clear maximum is found at InT ss 24.3(r rs 0.035THz). 
This is exactly the Zn(fii) at first minimum of Ct shown 
in Fig. 2. As shown before, this value corresponds to 
the typical frequency range occupied by one single res- 
onant peak. Such a result is coincident with our di- 
rect observation of the InT spectra that the most prob- 
able peaks/plateaus crossed by (InT) are those single 
peaks. Suppose the spectrum contains only sharp peaks 
but none plateaus. In(Acu) should be nearly Gaussian 
distributed with a single maximum at the mean value. 
However, the real probability distribution of InT shows 
a strange shoulder at InT « 26.6(r « OMITHz). Com- 



paring with Fig. 2 we find it exactly corresponds to the 
halfwidth of Ci n T, fis, which is just the typical width of 
SNS plateau, referring to the physical meaning of Ci 71 t- 
Actually, the measured probability distribution can be 
fitted very well by the summation of two Gaussian dis- 
tributions, as shown by the curves in Fig. 3. The pri- 
mary(left) Gaussian center gives the width of single res- 
onant peaks while the secondary (right) one gives exactly 
the typical width of the SNS plateau. Such a strange 
probability distribution implies substantial SNS plateaus 
with similar frequency widths contribute significantly to 
the probability distribution function of InT, resulting the 
characteristic shoulder of SNS. 

It is first time to see directly from the statistical distri- 
bution that the SNS is clearly distinguished from other 
localized states. From the distribution, we can see that 
the logarithmic plateau-width of SNS is Gaussian dis- 
tributed and its mean value agrees excellently with the 
results of correlation functions in third section, as ex- 
pected. So, both Ci n T and the probability distribution 
of InT provides proper physical values for characterizing 
the SNS. 



V. SYSTEM LENGTH DEPENDENCIES 

Both the correlation functions and the probability dis- 
tribution of InT suggest that the SNS favors a specific 
frequency width and a specific number (order) of local- 
ized states. But those studies are done with certain L/£. 
Next, to drive this conclusion further, we will calculate 
the correlation functions in the systems with different 
L/£. Fig. 4(a) shows Ct (the lower four blue curves) and 
CinT (the upper four red curves) functions for several L/£ 
values, where £ is fixed at 34.13a (W = 0.6) and L/£ = 
6, 10, 14, 18 from left to right. It clearly shows the peak 
width fix of resonant peaks of localized states (marked 
by blue arrows) decreases dramatically as increasing L 
while the decrease of plateau-width ils (halfwidth of Ci n T 
marked by red arrows) is much smaller, fix and ils ver- 
sus L/% are shown in Fig. 4(b) and (c). We can see fix 
decays exponentially with L (~ e~ L '£), as expected from 
traditional understanding^^. However, the decrease of 
ils is almost negligible and obviously much slower than 
fix- Actually, the plateau-width of SNS is determined 
by the intrinsic coupling strength(repulsing distance) be- 
tween localized states, which almost does not depend 
on L2*. Such a distinction naturally results a picture 
that, as increasing L more and more, the peaks of lo- 
calized states become so sharp that they are almost un- 
detectable, but the SNS plateaus are still there. In our 
numerical statistics, the number(order) of localized states 
increases nearly linearly with L/£. In Ref.— it is further 
argued that those SNSs (called intrinsic short necklace 
states) have main contribution to the fluctuation of trans- 
mission and also affect the value of localization length. 



VI. 



SUMMARY 



In summary, we have defined the basic quantities for 
SNS study and investigate the statistical properties of 
SNS in the strongly localized systems. We find two ap- 
proaches to quantitatively study SNS properties. The 
first approach is based on the correlation functions and 
the second one is based on direct measurements of the 
peak/plateau-width in logarithmic transmission spectra. 
In the first approach, we defined the correlation function 
CinT in InT spectra and show that the typical width of 
SNS plateaus can be characterized by the half-width of 
Ci n T- And, with the help of Ct (correlation function 
of transmission coefficient t), the most-probable order of 
SNS can be obtained. In the second approach, we de- 
fined the peak/plateau- width T in InT spectra and stud- 
ied the probability distribution of InT directly measured 
from a large number of spectra. The probability distri- 



bution of InT shows a novel shoulder and it can be fitted 
very well by the summation of two Gaussian distribu- 
tions. We showed that the first one is from the distribu- 
tion of localized state peaks and the second one is from 
the plateaus of SNS. The center of the second distribu- 
tion gives the average frequency width of SNS plateaus, 
which agrees very well with the value obtained from cor- 
relation function. As increasing the system length L, the 
plateau- width of SNS decays very slowly, compared with 
the exponentially-decayed peak- width of localized states. 
Finally, we note that the methods we used in our pa- 
per is not limited to the transport properties of localized 
system. They may be used for studying other statistical 
quantities which have similar properties with the trans- 
mission of localized systems. 
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